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Tiie optimum design, construction, and use of the LIGO, VIRGO, or LAGOS 
gravitational radiation detectors depends upon accurate calculations of their sensitiv- 
ity to different sources of radiation. Here I examine how to determine the sensitivity 
of these instruments to sources of gravitational radiation by considering the process 
by which data are analyzed in a noisy detector. The problem of detection (is a sig- 
nal present in the output of the detector?) is separated from that of measurement 
(what are the parameters that characterize the signal in the detector output?). By 
constructing the probability that the detector output is consistent with the presence 
of a signal, I show how to quantify the uncertainty that the output contains a signal 
and is not simply noise. Proceeding further, I construct the probability distribu- 
tion that the parameterization /i that characterizes the signal has a certain value. 
From the distribution and its mode I determine volumes V{P) in parameter space 
such that /2 G V{P) with probability P [owing to the random nature of the detector 
noise, the volumes V{P) are are always different, even for identical signals in the 
detector output], thus quantifying the uncertainty in the estimation of the signal 
parameterization. 

These techniques are suitable for analyzing the output of a noisy detector. If 
we are designing a detector, or determining the suitability of an existing detector 
for observing a new source, then we don't have detector output to analyze but are 
interested in the "most likely" response of the detector to a signal. I exploit the 
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techniques just described to determine the "most Ukely" volumes V{P) for detector 
output that would result in a parameter probability distribution with given mode. 

Finally, as an example, I apply these techniques to determine the anticipated 
sensitivity of the LIGO and LAGOS detectors to the gravitational radiation from a 
perturbed Kerr black hole. 



PACS numbers: 4.80.+z,98.60.a,97.60.Lf,6.20.Dk 
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I. INTRODUCTION 



Under the present schedule, both the United States Laser Interferometer Gravitational 
Wave Observatory (LIGO and the French/Italian VIRGO will begin operation 

in the late 1990s. Long before that time, theorists must lay a foundation for the study of 
gravitational radiation sources. Part of this foundation involves the construction of detailed, 
parameterized models of the waveforms from expected sources; another part involves the 
calculation of the anticipated sensitivity of the detector to each of these sources. Calculation 
of these kinds are not only needed for LIGO and VIRGO: design and technology studies 
for a Laser Gravitational-wave Observatory in Space (LAGOS) are currently being pursued 
@] and calculations of the sensitivity of LAGOS to appropriate sources are needed to guide 
these studies. 

In this paper I address the problem of calculating the anticipated sensitivity of a detector, 
like LIGO, VIRGO, or LAGOS, to an arbitrary source of gravitational radiation. The 
problem breaks up into two parts which I term detection and measurement. To "detect" is 
to decide whether the observed detector output contains a signal from a particular source 
or is just an example of noise; to "measure" is to assume the presence of a signal in the 
detector output and to characterize the signal in terms of the parameter(s) that describe 
the source (and its orientation with respect to the detector). 

Echeverria |P recently examined some of these issues in the particular context of deter- 
mining the precision with which one could characterize the mass and angular momentum 
of a perturbed Kerr black hole from observations in a gravitational radiation detector. The 
foundation of his analysis was the construction of a quantity similar to the signal-to-noise 
ratio (SNR), and he asserted that the parameters that characterize a signal observed in the 
output of the detector are those that maximized this quantity. This analysis is limited in 
two respects (Echeverria & Finn 0): 

1. The validity of the formalism is restricted to the limit of high SNR; and 
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2. The formalism cannot determine the amphtude of the signaL 

In addition, the conceptual basis of this calculation is not compelling: the determination of 
the parameters characterizing a signal in a noisy detector does not proceed by maximizing 
the SNR-like quantity defined by Echeverria [Q. 

In contrast, the techniques developed here are all based upon the construction of prob- 
abilities and probability densities. For the problem of detection, I construct the probability 
that the observed detector output is consistent with the presence (or absence) of a signal. 
In the case of the measurement problem, where the detector output is assumed to include a 
signal, the quantity constructed is the probability density that describes the hkelihood of a 
given signal parameterization. 

In §|I| I examine the twin processes of detection and measurement from the point of view 
of probability theory. The parameters characterizing a signal identified in the output of a 
noisy detector are defined to be those most likely to have resulted in the observed detector 
output. Some of the results described in this section are known elsewhere in the context of 
data analysis: they are included here for completeness and so that they may be compared 
with the techniques employed in Echeverria ||5[. In §|IT| I show how these same techniques 
can be exploited to evaluate the anticipated sensitivity of an instrument to a signal: i.e., 
how precisely can the parameterization of a signal observed in the detector be determined. 
I find both exact and, in the interesting limit of a strong signal, approximate techniques 
for evaluating the expected precision with which an observed signal can be described. As 
an example, in I apply the approximate techniques developed in §|T| to the determi- 
nation of the parameterization of the gravitational radiation from a perturbed black hole, 
especially the black hole mass M and dimensionless angular momentum parameter a. In §0 
I briefly compare the methods and results of Echeverria |Q with my own. My conclusions 
are presented in § [VI| 
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II. DETECTION, MEASUREMENT AND PROBABILITY 



In this section I consider two related problems that arise in the analysis of the output 
of a noisy detector: detection and measurement. The problem of detection is to determine 
whether or not a signal of known form [i.e., deterministic, though parameterized by one or 
several unknown parameters) is present in the detector output. The problem of measurement 
is to determine the values of some or all of the unknown parameters that characterize the 
observed signal. 

Note that the distinction between detection and measurement separates the determi- 
nation of the presence or absence of the signal from the determination of the parameters 
that characterize it: detection does not address the value of the unknown parameters, and 
measurement presumes the signal's presence. 

Detector noise can always conspire to appear as an example of the sought-for signal; 
alternatively, noise can mask the presence of a signal. In either case, noise interferes with 
our ability to determine the presence of the signal or the parameters that characterize it. 
Consequently, any claim of detection must be associated with a probability signifying the de- 
gree of certainty that the detected signal is not, in fact, an instance of noise. Similarly, when 
an observed signal is characterized it is appropriate to specify both a range of parameters 
and a probability that the signal parameters are in the given range. 

For example, I can examine the data stream from a gravitational radiation detector to 
determine (with some probability) whether the radiation from the I — \m\ — 2 mode of a 
perturbed, rotating black hole is present, irrespective of the black hole mass, angular mo- 
mentum, or orientation with respect to the detector. If I conclude that a signal is present 
in the data stream, then I can attempt to determine bounds on some or all of these param- 
eters, such that I expect the actual parameters characterizing the signal to fall within those 
bounds with a given probability. 

In the next several subsections I examine detection and measurement in more detail. I 
assume that the statistical properties of the detector noise are known, and also that the form 
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of the sought-for signal is known up to one or several parameters. My discussion focuses on 
determining the probability that a signal of known form is present in the output of a noisy 
detector, and on determining the probability that the unknown parameters have particular 
values. 

While the discussion in § [IV| is framed in the context of the measurement of gravitational 
radiation from astrophysical sources, the questions addressed in this (and the following) sec- 
tion are purely statistical ones and contain nothing that is specific to gravitational radiation, 
general relativity, or any particular physical system or theory. For more details, the reader 
may consult Wainstein & Zubakov 0. 



A. Detection 



Consider a data stream g{t) which represents the output of a detector. The data stream 
has a noise component n{t) and in addition may have a signal component m(t). The signal 
component is parameterized by several unknown parameters (denoted collectively as /x, and 
individually as /Xi); hence 



git) 



n{t) 



if signal not present. 



(2.1) 



n{t) + ■m{t; fi) if signal m(t; /x) present. 
Assume that /x is continuous, not discrete. I will describe how to determine the probability 
that m(t; fi), for undetermined /j,, is present in g{t), i.e., 



P{m\g) 



^The conditional probability that a^ 
signal of the form m{t; /x), for 
unknown fj,, is present given the 
^observed data stream g(t). 
Begin by using Baye's law of conditional probabilities to re-express P{m\g) as 

P{g\m)P{'m) 



(2.2) 



P{m\g) 



Pig) 



(2.3) 



where 



P{g\m) 



^The probability of measuring g ^ 



P(m) 



Pig) 



assuming the signal m is present. 



(2.4a) 



(2.4b) 



(2.4c) 



The a priori probability 
^that the signal m is present. 

The probability that the 

^data stream g{t) is observed. 

Also re-express P{g) in terms of the two possibilities m absent and m present, and further 
re-express the probability that m is present in terms of the probability that it is characterized 
by the particular fi: 

P{g) = P(^|0)P(0) + P{g\m)P{m) 

= P(y|0)P(0) + P(m) j dyp{tx)P[g\m{tx)], 

(2.5) 



where 



P(0)^ 



^(^10) 



P[g\m{i^)] 



^The a priori probability that ^ 



the signal is not present. 



^The probability density of ^ 



observing g[t) in the absence 
of the signal. 



p(m) = 



The probability density of 
observing g{t) assuming 
m{t] n) with particular 
is present. 

^The a priori probability ^ 

density that m{t) is 
y characterized by /i. 



(2.6a) 



(2.6b) 



(2.6c) 



(2.6d) 



Combining equations |2.3| and |2.5|, we find 



where 



^<"l«' = A + P(0)/P(,n) - 



A = I d^fxAifi) (2.8) 
M,).,i,)^. (2.0, 



In equation all of the dependence of P{m\g) on the data stream g has been gathered 
into the likelihood ratio A. Aside from A, P{m\g) depends only on the ratio of the a priori 
probabilities -P(O) and P{m). In turn, the likelihood ratio depends on two components: the 
a priori probability density p{m\fi) and the ratio P[g\m{iJ,)]/ P{g\0) . 

In order to determine P{m\g) we must assess the a priori probabilities and calculate 
the likelihood ratio. It is often the case that we know, or can make an educated guess 
regarding, the a priori probabilities. For example, the sources may be Poisson distributed 
in time [determining P(0) and P(m)], and they may be homogeneously distributed in space 
[determining p{r) oc r^, where r is the distance to the source]. At other times our assessment 
may be more subjective or based on imperfect knowledge, and in this case we can use the 
observed distribution of jj, to test the validity of our assessments using the techniques of 
hypothesis testing (Winkler [|| §7). 

Now turn to the evaluation of P[g\m{iJ,)]/ P{g\0). To determine this ratio, first note that 
the conditional probability of measuring g{t) if the particular signal m(t; n) is present is the 
same as the conditional probability of measuring g'{t) = g{t) — m{t; assuming that the 
signal m(t; n) is not present in g': 

P[g\m{fi)] = P[g-m{fi)\0]. (2.10) 

Consequently, we can focus on the conditional probability of measuring a data stream g{t) 
under the assumption that no signal is present [P(5'|0)]. 

In the absence of the signal, g{t) is simply an instance of n{t). Assume that the n{t) 
is a normal process with zero mean, characterized by the correlation function C„(r) [or, 
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equivalently, by the one-sided power spectral density (PSD) 5'„(/)]. In order to compute 
the ratio P[g\m{iJ,)]/ P{g\0) , consider the continuum hmit of the case of discretely sampled 
data {gi : z = 1, . . . , A^}, with the correspondence 



9i = giu), 

U - tj = it- j) At, 
T 



At 



A^- 1 



(2.11a) 
(2.11b) 
(2.11c) 



The probability that an individual g^ is a sampling of the random process n{t) is given by 

exp 



P{9i\0) 



1 9t 
2C„(0) 



[2vrC„(0)]^/' 

and the probability that the ordered set {g^ : i = 1, . . . N} is a sampling of n{t) is 



2.12 



^(^10) 



exp 



I]^fc=i Cjf}gjgk 



1 v^TV 
'2 



(27r)'^det||C„,, 



1/2 



(2.13) 



where Cj/^ is defined by 



and 



I 



(2.14) 



(2.15) 



(Mathews & Walker §14-6, Wainstein & Zubakov eqn. 31.11). Note that the normal- 
ization constant in the denominator of equation |2.13| is independent of the gf, consequently, 
it does not affect the ratio P[g\m{iJ,)]/ P{g\0) . Since it is this ratio that we are interested 
in, without loss of generality drop the normalization constant from the following. 
To evaluate equation |2.13| in the continuum limit, first note that 



6{tj-tk) = hm -^5jk. 



(2.16) 



Consequently, 



'jk 



1 /"OO /"OO 

lim — / dt^ e^-^*^ / dt, C^{t, - ti)C-\tu h 



At 



1 /"OO /"OO 

lim — / dtie^^'f'^C-\tutk) / rfre2-*^^a(r) 

At^O Al^ J-oo J-oo 

lim i^5„(/)C^i(/,4). 



(2.17a) 
(2.17b) 



To proceed from equation p.l7a| to p.l7b| , use the Wiener-Khintchine (c/. Kittel |jlO| §28) 
theorem to relate the PSD S'„(/) to the correlation function C„(r) and define 



2-Kift 



(2.18) 



Consequently, as we approach the continuum limit, we have 

n„2TTiftk 

C~^{f,tk) = lim Af 



(2.19) 



^*-o Snif) ■ 

With and Parseval's Theorem, we can evaluate the continuum limit of the argument of 
the exponential in equation |2.13| : 



N 



N 



l™o ^ Cjk9j9k= 1™ ^ ^^^^ \tj,tk)g{tj)g{tk 

T^oo j,k=l 



T^oo j,k=l 



lim -— 

Ai^O At 



At^O At2 



oo 



dtjdtkC ^itj,tk)gitj)g{tk) 
dfdhC^^{f,tk)nf)9{tk) 



-oo 
oo 



df 



Snilfl) 



(2.20) 



Sum 

Here and henceforth we will denote the Fourier transform of r{t) as r{f). 

Since the detector output g{t) is real, g*{f) = g{—f)- Define the symmetric inner product 

~9U)h\f) 



{g, h)^ df 



Sn{\f\) 



(2.21) 
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for real functions g and h. In terms of this inner product, 



Pi9\0) 

p{^l) exp [2 {g, m(/x)) - {m{^l),m{^l))] , (2.22) 



The hkehhood ratio A is found by substituting eqn. |2.22| into eqn 



To summarize, the probabihty P{m\g) that a signal of the class m{t; /x) is present in 
the output of the detector g(t) can be expressed in terms of three a priori probabilities 
[P(0), P{m), andp(/x)] and the ratio of two conditional probabilities [P{m\g) / P{0\g)]. The 
a priori probabilities must be assessed, while the ratio of the conditional probabilities can 
be calculated. Often we know, or can make an educated guess regarding, the a priori 
probabilities; at other times our assessment is subjective or otherwise based on imperfect 
knowledge. Finally we establish a threshold for P{m\g) [or, equivalently, for A, In A, or 
some other surrogate of P{m\g)], and say that if the P(m|5f) (or its surrogate) exceeds the 
threshold then we have detected the signal. 

I will not discuss detection further, except to say that the choice of threshold is influenced 
by our strategy to minimize errors. The two kinds of errors we can make are to claim the 
presence of a signal when one is in fact not present (a "false alarm"), or to dismiss an 
observed g{t) as noise when a signal is present (a "false dismissal"). In order to minimize 
the probability of a false alarm (conventionally denoted a) we want a large threshold, while 
to minimize the probability of a false dismissal (conventionally denoted f3) we want a small 
threshold. One obvious strategy for choosing the threshold is to minimize the sum a + /3, 
i.e., to minimize the probability of making an error. Alternatively, some other combination 
of a and f3 may be minimized, taking into account the relative seriousness of the different 
kinds of errors. Regardless, it is inadvisable to blindly choose a threshold for P{m\g) without 
careful consideration of the false alarm and false dismissal probabilities that arise and their 
relative severity. 
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B. Measurement 



Turn now to the question of measurement. From equations |2.3|, |2.5| and p]9| we have 



^The conditional probabihty that ^ 



p[m{fi)\g] 



the particular signal m{t] /j,) is 

present in the data stream g{t). 
A(/x) 



(2.23) 



A + P(0)/P(m) 

This conditional probability density is directly proportional to A(^t) and, since the denomi- 
nator in equation p.23| is independent of /x, it is maximized where A(^t) is maximized. If we 
assume that the signal is present, then the probability density that it is characterized by ^t 
is 



ml 



A 



(2.24) 



The goal of the measurement process is to determine a volume V{P) in parameter space 
such that n G V{P) with probability P. This volume is "centered" on the mode of the distri- 
bution p['m{fi)\g] in a way we define later on. The mode of either p[m{fi)\g] or p['m{fi)\g, m] 
is the n that maximizes A(/x). Denote the mode by /*.[] While I will occasionally refer to 
Ji as the "measured" parameterization of the signal, bear in mind that /i is only the most 
likely parameterization of the observed signal. 

If we assume that the global maximum of A(^t) is also a local extremum, then ft satisfies 

^A{^^) 







dfJ'i 



(2.25) 



equivalently, fi maximizes 



lnA(/x) = lnp(/x) + 2 {m{fi),g) - {m{^i) , m{^i)) 



(2.26) 



^While we assume in what follows that the distribution has a single mode, the generalization to 
a multi-modal distribution is trivial. 
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i.e., it satisfies 

This final set of equations is in general non-linear and may be satisfied by several different /x. 
Some will represent local maxima, while others will correspond to local minima or inflection 



points; thus, equation |2.25| is a necessary but not sufficient condition for p,. 

An important characterization of the strength of the signal in a detector is the signal- 
to- noise ratio (SNR). The "actual" SNR depends on the true parameterization of the signal 
Ji. We do not have access to /i; however, we do know that the most likely value of Ji is /i, 
and we define the SNR in terms of p,: 

p2 = 2(m(A),m(A)). (2.28) 

The factor of two arises because the power spectral density Sn{f) is one-sided while m(/) is 
two-sided. Note that is expressed in terms of the signal power {i.e., it is proportional to 
the square of the signal amplitude). There is some ambiguity in the literature over whether 
"SNR" refers to p or p^. We avoid the ambiguity by referring to either p or p^ wherever the 
context demands it. 

Having found the distribution p[m{^)\g\ (or p[m{^)\g ., m\) , we define the boundary of 
the volumes V{P) to be its iso-surfaces. The probability P corresponding to the iso-surface 
p[m(yLi)|(7, m] = K"^ is 

P = / d''pp[m{^l)\g]. (2.29) 

Note that since the distribution p[m(/x)|(yf, m] is not generally symmetric, p is not neces- 
sarily the mean of /x. Also, if the distribution p[m{^)\g\ has more than one local maximum 
then V{P) need not be simply connected. 

To summarize, suppose we have an observation g{t) which we assume (or conclude) 
includes a signal m{fi) (for unknown fi). We construct the probability density p[m{fi)\g,m\ 
according to equation |2.24| , and identify iso-surfaces of p[m{fi)\g,m] as the boundary of 
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probability volumes V{P) according to equation |2.29| . Finally, we assert that /x G V{P) 
with probability P. 



III. MEASUREMENT SENSITIVITY 

In §|I| we saw how to decide whether a signal is present or absent from the output of 
a noisy detector, and, if present, how to determine bounds on the parameterization of the 
signal. Now I show how to anticipate the precision with which a detector can place bounds 
on the parameterization that characterizes a signal. In particular, consider an observed 
g{t) which contains a signal m{t; Jx) for unknown Jx. We are interested ultimately in the 
distribution of 



S/j, = — fi 



(3.1) 



where jl is determined by the techniques discussed in §|I|. There are an infinity of possible 
g{t) that can lead to the same ft [corresponding to different instances of the noise n{t)], and 
for each there is a different probability distribution p[m{iJ,)\g] [cf. eqn. p.23| ) and a different 
set of probability volumes V{P). We will find the probability volumes V{P) corresponding 
to 



(3.2) 



'^The conditional probability density ^ 
that the signal parameterization is 
/i, assuming that the mode of the 
^distribution p[m{ii)\g] is fl. 

I show first how to do this exactly, and then show a useful approximation for strong signals. 
The mode p, of the distribution p[m{ii)\g^m\ satisfies 



(din \ 
(9 In n , _ , 



-2 i^in) 



(3.3) 
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;/. eqns. ^]T] and |2.27| ). Since n{t) is a normal variable with zero mean, so are each of the 



(n, dm/d^j) on the righthand side of equation |3]^. Denote these random variables Ui. 



n,3^(A)). (3.4) 

The joint distribution of the z/j is a multivariate Gaussian and its properties determine, 
through the equation the properties of the distribution of S/j,. Consequently we can 



focus on the joint distribution of the z/j. 

Since the z/j are normal, their distribution is determined completely by the means T7[, 
which vanish, and the quadratic moments 



To evaluate the average on the righthand side, we will use the ergodic principle to turn the 
ensemble average over the random process n into a time average over a particular instance of 
n. Recalling that a time translation affects the Fourier transform of a function by a change 
in phase, 

:p^r{t + T)]=e-^^'f^r[r{t)], (3.6) 



write 



(n(t + r),r(t)) = ^Jte-'-f^^^. (3.7) 



Consequently, 



1 /-^ 

(n, r) (n, s) = lim — / dr {n{t + r), r) {n{t + t),s) 



lim — / dr df y^ e-""^" / df ' V e 

T^oo2TJ-T J-oo ^ Snif) J-oo ^ S„{f') 

^f'-^ (3.8b) 

^{r,s). (3.8c) 
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In going from eqn. |3.8a| to eqn. |3.8b| , we used the definition of the PSD of the detector noise 

n{ty. 



hm i|n(/)|^ 

1 — >oo ± 



(3.9) 



cf. Kittel [10] §28). With the result in equation [3.8c|, we have 



/ din \ / diTi 



I dm dm 



= a 



-1 



(3.10a) 
(3.10b) 



In terms of the Cij {i.e., the inverse of the Q^^), the joint distribution of the z/j is given by 



exp 



2 Si J CijViVj 



1/2 



(3.11) 



{2Tir detllQ, 

This is also the joint distribution of the quantities 

-2(,n(M)-,n(A).|^(w)-^(M) 

that appear on the lefthand side of equation |3]^; consequently, we expect that for an obser- 
vation characterized by a given p, the probability volumes V{P) are given implicitly by 

d\np 



(3.12) 



X 



2{m{fi) -m{fi),—{fi) 

(din 
m{fi) -m{li),—{fl) 



dlnp 



(3.13) 



where 



exp 



2 X^i.j CijPiUj 



Y...,<^^,^^.<^' [(27r)^det||Qi 



1/2- 



(3.14) 



This result is exact as long as the maximum /i of A(^) is also a local extremum of A(^). 

As the SNR becomes large the distribution p(/i|/i) becomes sharply peaked about /i and 
the determination of the volume V^(-P) is greatly simplified. Suppose that (? is so large that 
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for fi G V{P) for all P of interest, the difference m(/j,) — m{fl) can be linearized in Sfi. We 



then obtain in place of equation 3.3 



The random variables 5// are related to the by a linear transformation, 



(3.15) 



5 Hi 



d\n.p 



(3.16) 



consequently, the S^i are normal with means 



and quadratic moments 



^ (9 In J} 



(3.17) 



(3.18) 



The probability distribution p{6iJ,\fi) is a multivariate Gaussian (c/. eqn. |3.11| ): 



exp 



{Sf^i - Sf^i) {Sf^j - 



(27r)^det||Ci_ 



1/2 



(3.19) 



Note that the matrix Cij now has acquired a physical meaning: in particular, we see that 
the variances of the 5 Hi are 



erf = (6Hi - 8lJL, 

= C 

and the correlation coefficients rj,- are given by 



(3.20) 



C 



(3.21) 



In this sense we say that Cij is the covariance matrix of the random variables 6fi. 

In the strong signal approximation, the surfaces bounding the volume V{P) are ellipsoids 
defined by the equation 
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(3.22) 



i,3 



where the constant K is related to P by 



exp 



2 ^—ii,] ^ij 



(27r)^det||Ci_ 



1/2 



(3.23) 



It is often the case that not all of the parameters that characterize the signal are of 
physical interest. In that case, we may integrate the probability distribution (eqn. p. 11 



or 



|3.19|) over the uninteresting parameters, leaving a distribution describing just the parameters 
of physical interest. 



Finally we come to the question of when the linearization in equation |3.15| is a reasonable 
approximation. Two considerations enter here: 

1. It is important that the probability contours of interest {e.g., 90%) do not involve Sfi 
so large that the linearization of m{Ji) — m{ji) is a poor approximation; and 



2. It is important that the condition number (c/. Golub & Van Loan |TT[]) of the matrix 
Cj^^ be sufficiently small that the inverse Cij is insensitive to this approximation in the 
neighborhood of /i.0 

These two conditions will depend on the problem addressed. If the validity of the lineariza- 
tion procedure is doubtful owing to the violation of either or both of these conditions, then 



we must fall-back on equation 3.3 and the exact results in equations 3.13 and 3.14 



IV. APPLICATION: A PERTURBED BLACK HOLE 

In this section, I show how to use the approximate techniques developed in § [m| to find 
the precision with which the mass and angular momentum of a perturbed black hole can be 

^Recall that the relative error in is the condition number times the relative error in C^j^: for 
a large condition number, small errors in C^-^ introduced by the linearization approximation can 
result in large errors in 5^i. 
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determined through measurement in an interferometric gravitational wave detector. This 
problem was first considered by Echeverria 0. 

Consider a single interferometric gravitational wave detector and a perturbed black hole 
of mass M and dimensionless angular momentum parameter a. Focus attention on a single 
oscillation mode of the black hole, e.g., the I = m = 2 mode. The strain measured by the 
detector has the time dependence of an exponentially damped sinusoid characterized by the 
four parameters Q, f, V, and T: 

for t < 0, 

h{t) = (4.1) 
1^ y^i/3^-.fit-T)/Q gij^ [27r/(t - T)] for t > 0. 

For convenience, assume that the perturbation begins abruptly at the starting time T. The 
frequency f depends inversely on the mass of the black hole, and has a weak dependence on 
its angular momentum: for the I = m = 2 quasi-normal oscillation mode, 

/ - ^ (4.2a) 
F{a)^l-—{l-af'\ (4.2b) 

is an accurate semi-empirical expression for the real part of the quasi-normal mode frequency 
(Echeverria eqn. 4.4 and tbl. II). The quality Q is the damping time r measured in units 
of the frequency /: 

Q = nfr. (4.3) 

For the I = m = 2 oscillation mode of the black hole, Q depends entirely on a: 

Q ~ Q(a) = 2(1 - a)-9/20 (4.4) 

(Echeverria eqn. 4.3 and tbl. II). Finally, the amplitude V^^^'^ of the waveform depends 
on the distance to the source, the size of the perturbation, and the relative orientation of 
the detector and the source. 

This peculiar parameterization of the amplitude refiects our expectation that perturbed 
black holes are distributed uniformly throughout space {i.e., V oc r^) and that all relative 
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orientations of the detector and the black hole source are equally probable. Additionally, 
it reflects an assumption that perturbations of any allowed amplitude are equally probable; 
consequently, the a priori distribution piV) is uniform. Let us also assume that p(a), p{f) 
and p{T) are uniform and that there is no a priori correlation of a, /, V , or T. 

An interferometric gravitational wave detector is naturally a broadband receiver, though 



it can be operated in a narrow band mode (c/. Vinet, Meers, Man, & Brillet [0, Meers |T3 



and Krolak, Lobo & Meers IQ)- Assume that the detector response function is uniform in 



the frequency domain over the bandwidth of the gravitational wave; consequently, the signal 
component in the output of the detector [m(t, ^)] is equal to the waveform h{t] Q, /, V, T) {cf. 
eqn. [4.1|) . Assume also that the noise PSD (Sn) of the detector is independent of frequency 
(/) in the bandwidth (1/r) of the signal (I will discuss the validity of this approximation 
below) . 

A. The signal-to-noise ratio 

As a first step toward finding the precision with which a, M, V and T can be determined. 



we calculate the SNR p^. With h given by equation ^Aj evaluate using equation |2.28| to 
obtain 

^ " 7rfVy^{l + AQ^)Sn ^^'^^ 

This expression is valid to better than a percent as long as the signal is observed for a period 
of time At > 2.5r. 

In arriving at equation [4. 5| we assumed that the noise PSD is constant over the bandwidth 
of the signal so that Sn = Sh{f)- The signal bandwidth Af is approximately 

A/-- = ^(l-af/^V. (4.6) 

For small a the bandwidth is approximately /, while for large a the signal is monochromatic. 
For small a the approximation that Sn is constant over the bandwidth of the signal is only 
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a fair approximation for LIGO (c/. Vogt [|l|, Abramovici et al. p[) or LAGOS (c/. Faller 
et al. 0]); however, it becomes a good approximation for both detectors when a > 0.9 
(corresponding to Af / f < 1/2). 

The amphtude V depends on the detector orientation with respect to the black hole, the 
amplitude of the perturbation, and the distance between the black hole and the detector. 
Average over all possible orientations of the detector with respect to the black hole (c/. 



Thorne |jT5| §9.5.3) to obtain 



^ 16 

p2 - — 



eM /My 



(4.7) 



5 F2 (1 + 4Q2) Sn \ r J 
where eM is the total energy radiated by the I = m = 2 mode of the black hole perturbation 
and r is the distance of the source. 

When operated as a broadband detector, the LIGO advanced detectors will be most 
sensitive to perturbed black holes with 50 ^ M < 100 where Sn — 10~^'^Hz~^ (c/. 
Krolak, Lobo & Meers fl^, Dhurandhar, Krolak & Lobo |jl6|, Vogt et al. Abramovici et 
al. 0). LAGOS will be most sensitive to perturbed black holes in the range 10^ Mq ^ M < 
WMq, where ~ 10^^^ Rz'^ (cf. Faller et al. [|]). Consequently 



—1/2 
p2 



5.8^(0) < 



( e / 3Mpc \ / M Y'^ f ] 

V4X10-V \ r J VsOMo^ V S„ 



V 5x10-3 J \ r J VlO'^MeJ V" 



LIGO 
LAGOS 



(4.8) 



where 



G{a) = 



37 
200 



17Q2 



1/2 



(4.9) 



F2(l + 4g2) 

For frequencies outside of the range 100 Hz to 200 Hz, the LIGO PSD Sn scales with fre- 
quency: for frequencies greater than approximately 200 Hz (corresponding to M < 50 Mq), 
Sn scales as (c/. Thorne [|I^, Krolak, Lobo & Meers Jl^)) and for frequencies less than 
100 Hz (M > 100 Mq) it scales as /"^ (c/. Dhurandhar, Krolak, & Lobo |g). Similarly, for 
frequencies outside the range 10"^ Hz to 10^^ Hz the LAGOS PSD Sn scales with frequency: 
for / > 10"^ Hz it scales as p, and for / < 10"^ Hz it scales as /^'^ {of. Faller et al. g]). 
Consequently 
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_ M>100Mo 

p2 oc { (4.10) 

M < 50 M0 



for LIGO and 



_ M-i M>1O^M0 

p2 oc <{ (4.11) 

M < 10^ Mq 

for LAGOS. 

Little is known about the rate of, or the energy radiated during, black hole formation [cf. 
Rees |]T^, Kochanek, Shapiro & Teukolsky [|18|); however, owing to the extreme sensitivity 



of both the LAGOS and LIGO detectors, it seems a conservative estimate that the formation 
of a black hole of mass 10^ Mq ^ M ^ 10"^ Mq anywhere in the universe will be detectable 
by LAGOS, and the formation of black holes with 50 Mq ^ M ^ lOOM© will be observable 
in LIGO at least to the distance of the Virgo cluster (~10 Mpc). Additionally, note that 
the energy radiated in the / = 2 mode during the radial infall of a test body (mass m) onto 
a Schwarzschild black hole (mass M) is given by 

2 

771 

= eM ~ IQ-^ — , (4.12) 



(Davis, Ruffini, Press, & Price |19|, Oohara & Nakamura |20|; similar results hold for Kerr 



black holes: Sasaki & Nakamura Kojima & Nakamura |^). Consequently, the capture 
of a solar mass compact object {e.g., a black hole or neutron star) onto a black hole of 
mass 10^-10'' Mq (corresponding to e ~ 10"^^-10^^^) may also be observable to a distance 
of 3 Mpc {cf. eqn. ^) 

Figure |I| shows the factor G{a) {cf. eqn. ^4.9| ) as a function of a. This figure may also be 
regarded as a plot of p{a) for fixed M, e, r, Sn, and detector-source orientation. With this 
interpretation, note how p decreases with increasing a. The reason for this behavior is that 
at fixed M, the frequency / and damping timescale r both increase with a; consequently, a 
signal of smaller amplitude {i.e., smaller p^) will yield the same radiated energy. 
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B. Precision of measurement 



While the parameters Q and / are convenient for characterizing the detector response, 
it is the determination of a and M that is of direct physical interest. If the perturbed 
black hole is also observed electromagnetically {e.g., if it is the result of the gravitational 
collapse of a star in a type II supernova), then determination of V and T may also be 
interesting. Regardless, we are more interested in the covariance matrix for the parameters 
{a, M, V, T} than for the parameters {Q, /, V, T}. It turns out, however, that it is simpler 
to first determine the covariance matrix for the parameterization {Q, /, V, T}. 

To find the covariance matrix for the parameterization {a, M, V, T}, first define the 
three dimensionless parameters e', ^' and C by 



(4.13a) 



Vi' = V- V, 



(4.13b) 




(4.13c) 



and evaluate C[, ^ for the parameterization {Q, e', C,'}: 




(4.14) 



The components of C- ^ appearing in equation |4.14| are 




(4.15d) 



(4.15b) 



(4.15a) 



(4.15c) 
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The covariance matrix Cij for the parameterization {a, e, ^, (} is given in terms of C-j by 



k,l 



(4.18) 



where the symmetric matrix J'ij is given by 



J 



\ 



dQ n n 

da fM da ^ " 

{fMf 

1 

ifMY 



(4.19) 



Like C^j, the matrix Cij is a function only of a and p^, and has the elements 



^ae 
Ca( 



(l + 2Q^)(l + 4g^) 1 

/ [QF^ (1 + 2g^) (1 + 4Q^) - 2FQ'] F' 

1-2Q^ + 8Q^ \ 1_ 
+ 2g4(l + 4g2) j p2' 

9(i + 2g^)^ 
2gv ' 

8Q' 



(4.20a) 



C. 



(i + 4g2)FV' 

g (1 + 2Q^) (1 + 4g2) - Fg^ 

2Fg3g'V2 
3(i + 4g2)(i + g2) 

2g3g'V2 

3 [g (1 + g^) (1 + 4g2) + (g^ - 1) Fg^] 
2g4Fgy 

(i-4g2)Fg^-g(i + 4g2)F^ 
g(i + 4g2)FW 



gFp2 



(4.20b) 

(4.20c) 

(4.20d) 

(4.20e) 

(4.20f) 
(4.20g) 

(4.20h) 
(4.20i) 

(4.20j) 



where F(a) and Q{a) are given by equations |4.2b| and [4.4| . Finally, in terms of these coeffi- 
cients, we have (c/. eqns. |3.20| and p.21| ) 
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2 71 /r2 2 

aM = M 


(4.21a) 




/a r\ -1 1 \ 

(4.21b) 


2 71 /r2 2 

erf, = M (T^ 


(4.21c) 


'"aM = fae 


(4.21d) 


TaV = 


(4.21e) 


TaT = Tad 


(4.21i) 


tmv = 


(4.21g) 




(4.21h) 


rvT = r-^c- 


(4.21i) 



The results for the standard deviations cr^ and ctm/M and correlation coefficient 
found semi-numerically in Echeverria |P (his eqns. 4.10a-c and table II) are approximations 
to the analytic results found here in equations 4.20a , 4.21a and 4.21d| . Additionally, we give 
analytic forms of the other variances and correlations coefficients. 

Figure shows (Ta, ctm/M, ay /V , and ctt/M (eqns. [4.2Ua| - ^I^U^ ), normalized by 
p as shown, as functions of the measured a. Note how for a ^ 0.8, the angular momentum 
parameter is determined less precisely than the mass. Figures |]a-c shows the six corre- 
lation coefficients (c/. eqns. p.21| , [4.20eH4720]l) as functions of the measured a. These are 
independent of p. In figures ^-c note how 5a and 5M are highly correlated so and are not 
statistically independent parameters: for a complete discussion of this point, see Echeverria 

i§iv. 

Figure ^ shows the standard deviations for fixed p. It is also useful to consider these same 
quantities for fixed M, r, e and Sn as was done in equation p7| , |47^ and figure |l|. Defining 
Po by 



p = poG{a) 



(4.22) 



where G{a) is given in equation figure ^ shows cTq, aM/M, oy jV and ax/M normalized 
by Po (for LAGOS and LIGO, po is given by equation [4.8|) and as functions of a. Note the 



26 



difference between figures |] and ^: in the first case, the SNR is held constant while in the 
second the the energy of the perturbation is held constant. In the second case, the precision 
with which a and M can be determined does not increase as rapidly with a as in the first 
case, the precision with which T can be determined is independent of a, and the precision 
with which the amplitude can be determined decreases with increasing a. 

The elements of Cj^^ fully determine the distribution p[m[^)\g\ and the volumes V{P). 
Generally we will have no interest in V and T, in which case we integrate the distribution 
over all T and V to find the two-dimensional distribution 

p[m{a,M)\g\ = 



exp 



2 1- 



AM I Aa 



) AaAMrg 

O-aO-M 



aM) 



where 



(4.23) 



AM = M-M 



Aa = a — a. 



(4.24a) 
(4.24b) 



V. DISCUSSION 

The earlier results of Echeverria on the precision of measurement are restricted to the 
case of large p where the distribution of 6fi is well-approximated by a Gaussian, though there 
is no discussion of what constitutes a sufficiently large SNR. Additionally, those results do 
not provide any guidance for estimating the precision with which the amplitude of the signal 
can be measured. Finally, there is no clear connection drawn between the measurement of 
the fl, the estimates of 6fii6fij, and the probability that |/i — /ip < 6pf. The restriction 
to strong signals is required because of the expansion of Echeverria's expression for p in 
a power-series about the "measured" parameters and also because the methods described 
fail to take into account prior knowledge about the distribution of the parameters [i.e., the 
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This prior knowledge plays an important role when the distribution of the S^i is not 
uniform or the SNR is small. 

On the other hand, the maximum likelihood analysis described in §^ is applicable for all 
SNR (though the approximate techniques discussed at the end of § PTT| are appropriate only 
when the distribution of the 5^ is well-approximated by a Gaussian). It does not elevate any 
parameter to a special status: the amplitude of the signal and its precision are determined 
in the same way that all other signal parameters and their precision are determined. Finally, 
it gives clear meaning to the measured parameters p, and the precision of measurement by 
providing the probability distribution of the 5^1. 

A complete discussion of how the optimal filter techniques of Echeverria are related 
to the maximum likelihood analysis presented here can be found in Echeverria & Finn 0. 

VI. CONCLUSIONS 

In the analysis of the results of an observation made in a detector, it is useful to distin- 
guish between detection and measurement. The analysis involved in detection refers only 
to the presence or absence of a signal characterizing a particular class of sources to which 
the detector is sensitive {e.g., perturbed black holes). A particular source of this class is 
described by a set of parameters: e.g., among the parameters describing the signal from a 
perturbed black hole is the black hole mass and angular momentum. Detection addresses 
only whether a signal of this class is present in the observed output of the detector, and not 
the particular values of the parameters that best describe the signal. 

Measurement follows detection: it refers to the determination of the values of the pa- 
rameters that best characterize the particular signal assumed to he present in the detector 
output (it only makes sense to speak of measuring the parameters of a real signal). For 
example, once we have concluded that we have detected the signal from the formation of a 
black hole, then we can go on to measure the black hole mass and angular momentum. 

In order to determine whether the observed output of a detector includes a signal from 
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a given class of sources, we saw how to calculate the probability that the detector output is 
consistent with the presence of the signal. That probability depends on the characteristics of 
the detector noise, the observed detector output, and a parameterized model of the detector 
response to the signal. In addition, it depends on several a priori probabilities that must 
be assessed. When the calculated probability exceeds a certain threshold then we say that 
the we have detected a signal. Setting the threshold requires careful consideration of the 
relative severity of falsely claiming a detection and incorrectly rejecting a signal. 

To determine the values of the parameters that characterize the detected signal, we saw 
how to construct the probability distribution that describes how likely different parameter- 
izations /x are. We identified p, as the mode of the distribution, i.e., the parameterization 
that maximized the probability density, or the most likely parameterization. Owing to de- 
tector noise, p, differs in a random fashion from the unknown Ji that actually describes the 
signal. We characterized our uncertainty over the actual description of the signal by speci- 
fying a volume V{P) in the parameter phase space, centered on /i, such that the Ji G V{P) 
with probability P. 

We then proceeded to exploit these techniques to anticipate the precision with which 
the parameterization of a particular signal can be determined by a given detector: i.e., we 
evaluated the sensitivity of the detector to the signal from a class of sources. 

To do so, we found the probability distribution of ft — jl and defined volumes V{P) in 
phase space such that fi G V{P) with probability P. These volumes determine the precision 
with which we expect we can determine the signal parameters in a real observation. In 
the interesting limit of a strong signal the anticipated probability distribution oi Ji — p, for 
fixed p is close to Gaussian and the associated volumes V{P) are ellipsoids. In this limit 
we found approximate techniques for determining the size and orientation of this ellipsoid. 
Both the exact and approximate expressions developed provide a powerful means of studying 
the sensitivity of a proposed detector or detector configuration to a source of gravitational 
radiation. These techniques are currently being employed to study the sensitivity of the 
LIGO detectors to binary coalescence precessing axisymmetric neutron stars and 
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non-axisymmetric neutron stars p5 



As an example of the process of measurement, we evaluated the variance in the mass 
and angular momentum of a perturbed black hole as determined by observations in a gravi- 
tational wave detector. These results improve upon those found earlier [cf. Echeverria P), 
and we discussed the origin of the differences. 

The LIGO detector, currently under construction, and the LAGOS detector, currently 
being designed, are both very sensitive to gravitational radiation from perturbed black holes. 
A perturbation of a 50 Mq-IOO Mq black hole that radiates as little as 10~^ of the black hole 
mass should be observable with LIGO at the distance of the Virgo cluster of galaxies, and a 
perturbation of a 10^ Mq-IO"^ Mq black hole that radiates as little as 10~^ of the black hole 
mass should be observable by LAGOS throughout the Universe {cf. eqn. ^I8|) . 
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FIGURES 

FIG. 1. The expected signal-to- noise ratio (SNR) of the I = m = 2 mode of a perturbed black 
hole as a function of the angular momentum parameter a. The dependence of the SNR on the 
black hole mass, distance, total energy radiated, and the detector noise PSD has been scaled out, 
leaving only the dependence on the angular momentum parameter. For more details see equation 



4.8 and surrounding text. 



FIG. 2. The expected standard deviation of the black hole angular momentum parameter (da), 
mass (aM/M), initial moment of perturbation [gt/M), and perturbation amplitude {ay /V) as a 
function of the angular momentum parameter a. The dependence of these standard deviations on 
the signal-to-noise ratio (SNR) p has been scaled out as shown. For more details see equations 



4.20a, 4.21a-4.21c and surrounding text. 



FIG. 3. The correlation coefficients for errors in the angular momentum parameter a, mass 
M, initial moment of perturbation T, and perturbation amplitude y as a function of angular 
momentum parameter. For more details see equations [4.21d - 4".21i| and surrounding text. 



FIG. 4. Like figure 2, except that the total energy radiated by the perturbation is held fixed 
instead of the SNR. Compare with figures |^ and |2|. For more details see |IV[ . 
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